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We present a detailed study of the quark and lepton mass spectra in a 50(10) framework with 
one IQh and one 12Qh Higgs representations in the Yukawa sector. We consider in full generality 
the interplay between type-I and type-II seesaw for neutrino masses. We first perform a \ 2 fit of 
fermion masses independent on the detailed structure of the GUT Higgs potential and determine 
the regions of the parameter space that are preferred by the fermion mass sum rules. We then apply 
our study to the case of the minimal renormalizable SUSY 5*0(10) GUT with one 10.tr, one 126h, 
one 126//, and one 210 h Higgs representations. By requiring that proton decay bounds are fulfilled 
we identify a very limited area in the parameter space where all fermion data are consistently 
reproduced. We find that in all cases gauge coupling unification in the supersymmetric scenario 
turns out to be severely affected by the presence of lighter than GUT (albeit B — L conserving) 
states. We then conclusively show that the minimal supersymmetric S'O(IO) scenario here considered 
is not consistent with data. The fit of neutrino masses with type-I and type-II seesaws within a 
renormalizable 50(10) framework strongly suggests a non-SUSY scenario for gauge unification. 
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I. INTRODUCTION 



Understanding the pattern of fermion masses and 
mixings is one of the longstanding challenges in par- 
ticle physics. In this respect Grand Unified The- 
ories (GUTs) do provide an appealing and power- 
ful tool to address the multiplicity of matter states, 
by involving stringent relations among the known 
fermions (and possibly implying a natural enlarge- 
ment of the minimal sector). Very appealing candi- 
dates for a GUT are models based on the SO(10) 
gauge group All the known fermions plus right- 
handed neutrinos fit into three 16-dimensional spino- 
rial representations of 5*0(10), and hence the model 
naturally leads to neutrino masses based on the see- 
saw mechanism [H, 0] ■ 

In this work we concentrate on a minimal renor- 
malizable version of the supersymmetric SO(10) 
GUT, where fermion mass matrices are obtained 
from the Yukawa couplings of the matter fields to 
10 and 126 dimensional Higgs representations Q. 
The 126// representation contains the scalar multi- 
plets (10, 3, 1) © (10, 1, 3) under the Pati-Salam sub- 
group 1 SU(4)c®SU{2) L ®SU(2) R , whose vacuum 
expectation values (VEVs) generate mass matrices 
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for the left- and right-handed neutrino fields respec- 
tively. Hence, in general there are contributions to 
the effective mass matrix of the light neutrinos from 
the type-I as well as type-II [3| seesaw mecha- 
nisms. Thanks to its minimality the model leads to 
constrained relations among quark and lepton mass 
matrices @, 0|: 

M d = vfY 10 + vf 6 Y 126 , 

M u = v]?Y w + v™Y 126 , 

Mi = vfY w - 3vf 6 Y 126 , (1) 

M D = vl°Y w - 3v™Y 126 , 

M L = v L Y 126 , 

Mr = V R Yi 26 . 

Here Mr>, Mr,, Mu denote the Dirac neutrino mass 
matrix, the mass matrix of the light left-handed 
neutrinos, and the mass matrix of the heavy right- 
handed neutrinos respectively. The indices d, u refer 
to down- and up-quarks, I denotes the charged lep- 
tons, while Yiq and Y\ 2 q are two complex symmetric 
matrices related to the 10// and 126// Yukawa inter- 
actions. The various v's denote the VEVs of the rel- 
evant Higgs multiplets, with v^ d 126 standing for the 
SU(2)l doublet components giving rise to the light 
MSSM-like Higgses, while v R (v L ) are the SU(2) L 
singlet (triplet) VEVs entering the type-I (type-II) 
seesaw formulae for the light neutrino masses. The 
effective mass matrix for the light neutrinos is then 
written as 

M v = M L - MdM^Md , (2) 
where the symmetry of Mrj is used. In Eq. flU the 
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first and second terms correspond to type-II and 
type-I seesaws respectively. Hence, all fermion mass 
matrices are predicted in terms of two complex sym- 
metric Yukawa matrices and six VEVs. 

Quite an amount of activity has been devoted re- 
cently to the study of the minimal renormalizable 
SUSY SO (10) model, see for instance Refs. |M [9, 
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231 \24 . [2a . [26j . This interest has been triggered to 
some extent by the experimental progress in neu- 
trino physics, and the observation of Refs. [H ITlj . 
that within the minimal 50(10) model with type- 
II seesaw dominance b — r unification allows natu- 
rally for a large lepton mixing due to the destructive 
interference between the 33 elements of the down- 
quark and charged lepton mass matrices. Analy- 
ses of proton decay within this setting are found, 
e.g., in Refs. [13, S|, while various extensions of 
the model have been considered in view of rising 
some of the tensions with the data, by adding non- 
renormalizable terms [l7| . by minimally extending 
the Higgs sector [l8|, or by including Yukawa cou- 
plings of the fermion fields with an additional 120- 
dimensional Higgs representation [H [24|, [2(| . 

Most of the previously quoted analyses consider 
the dominance of one type of seesaw and assume 
that the vacuum gives the correct neutrino mass 
scale. Very recently two studies appeared show- 
ing an intrinsic antagonism between seesaw neu- 
trino masses and coupling unification in the minimal 
SUSY 50(10) scenario [H[3(|, pointing out a criti- 
cal tension between the needed seesaw neutrino scale 
and the detailed spectrum arising from the minimal 
50(10) breaking GUT vacuum UMIH, [H|]. 

The aim of the present paper is twofold. First, 
in Sees. |TT] and IIIII we study in detail the implica- 
tions of the fermion mass sum rules that emerge 
from 10 H and 126 H Yukawa terms in a 50(10) GUT 
framework. The analysis is performed in full gener- 
ality, allowing for complex Yukawas and VEVs, and 
considering the neutrino mass matrix as originat- 
ing from an admixture of type-I and type-II seesaw 
as given by Eq. In the second part of the pa- 
per, Sec. IIV[ we specialize our study to the minimal 
renormalizable SUSY 5O(10) model, by including 
the constraints coming from the detailed structure of 
the GUT symmetry breaking vacuum and perform- 
ing an exhaustive study of the parameter space by 
optimizing a % 2 function. This approach is comple- 
mentary to the random parameter searches applied 
in previous studies, since the algorithm converges to 
an optimal solution (if it exists), and isolated solu- 
tions (allowed by an acceptable fine tuning of the 
parameters) can be found, which are easily missed 
in a random parameter scan. 

We show conclusively that the minimal renormal- 
izable 5O(10) supersymmetric scenario does not al- 
low for a consistent fit of the fermion mass spectrum, 



the model being unable to reproduce the correct neu- 
trino mass scale via type-I and/or type-II seesaw. 
Our negative result can be cast in more general (and 
simple) terms considering that, while i) the present 
constraints on proton decay force any GUT scale 
to lie above 10 16 GeV and ii) SUSY gauge unifica- 
tion works very well considering a desert scenario 
between the weak and the GUT scales, the neu- 
trino mass coming from type-I and/or type-II see- 
saws, being proportional via Yukawa and Higgs po- 
tential couplings to w^, eafc /M, requires for m wea f. ss 
200 GeV, M < 10 15 GeV. As a consequence (barring 
a strong interacting sector in the theory), interme- 
diate mass scales appear which are bound to affect 
SUSY gauge coupling unification at the least. 

All this applies to the minimal renormalizable 
setup. As it was pointed out very recently [32], non- 
minimal realizations of the renormalizable Yukawa 
sector such as those containing an additional 120- 
dimensional Higgs multiplet may provide a way 
out of the issues: the neutrino mass scale can be 
enhanced by the type-I seesaw with right-handed 
neutrino masses below the GUT scale due to tiny 
Yukawa couplings associated to the 126# multiplet, 
whose role in the charged matter sector can be 
taken over by the new 120jy multiplet. Although 
the simplest realization of this programme (where 
the charged matter sector is dominated by 10# and 
120 h and the Yukawas of 126# are neglected) seems 
to fail [HI , this attempt does provide a direction for 
further studies. Alternatively, Planck induced non- 
renormalizable operators may also provide the scale 
suppression needed by the neutrino sector, see for 
instance Ref. [34[ and references therein. 

On the other hand, the simplicity and the high 
level of predictivity that makes the minimal renor- 
malizable SUSY 5O(10) GUT extremely appealing 
is lost in both extensions. Before calling the minimal 
setup a dead end an extensive analysis of potential 
loopholes in the abovementioned arguments is due. 
This aim we pursue with the present paper up to, 
helas, the bitter end. 



II. DESCRIPTION OF THE ANALYSIS 

A. The parameterization 

In order to investigate whether Eqs. (JTJ) and ([2]) al- 
low for fermion masses and mixing in agreement with 
the data we proceed as follows. It turns out to be 
convenient to express the Yio and Yn§ Yukawa ma- 
trices in terms of Mi and M d , and substitute them 
in the expressions for M u , Mo and M v : 

M u = f u [(3 + r)M d + (1 - r)M e ] , (3) 
M D = f u [3(1 - r)M d + (1 + 3r)M e ] , (4) 
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where 



observable 



input data 



fu 



1 v. 



10 



4 V 10 > <0 W 126 ■ 

The neutrino mass matrix is obtained as 



(5) 



M v = f„ 



(M d - Mi) 



,M D , , , Mjj 

Ju JU 



with 



_l_vr_ 

J» 4 v 126 



(6) 



, (jUt , (7) 



The parameter |£| controls the relative importance 
of the type-I and type-II seesaw terms: For |£| — > 
one obtains pure type-II seesaw, whereas |^| — s- oo 
(with — const.) corresponds to type-I seesaw. 
For |£| ~ 1 both contributions are comparable. 

In what follows we denote diagonal mass matrices 
by m x , x = u,d, i, v, with eigenvalues corresponding 
to the particle masses, i.e., being real and positive. 
We choose a basis where the down-quark matrix is 
diagonal: Md = rhd- In this basis Me is a general 
complex symmetric matrix, that can be written as 
Mi = WjfhfWg, where We is a general unitary ma- 
trix. Without loss of generality /„ and /„ can be 
taken to be real and positive. Hence, the indepen- 
dent parameters are given by 3 down-quark masses, 
3 charged lepton masses, 3 angles and 6 phases in 
fu, fu, together with two complex parameters 
r and £: 21 real parameters in total, among which 
8 phases. Using Eqs. ([3]), (fj|, and © all observ- 
ables (6 quark masses, 3 CKM angles, 1 CKM phase, 
3 charged lepton masses, 2 neutrino mass-squared 
differences, the mass of the lightest neutrino, and 
3 PMNS angles, 19 quantities altogether) can be 
calculated in terms of these input parameters. Al- 
though the number of parameters is larger than the 
number of observabels the system is sensibly over- 
constrained due to the non-linear structure of the 
problem. 

Since we work in a basis where the down-quark 
mass matrix is diagonal the CKM matrix is given by 
the unitary matrix diagonalizing the up-quark mass 
matrix up to diagonal phase matrices: 



W U M U W* 



(8) 



with 



W u =diag(e ipi ,e 



)^ C KMdiag(e Ml ,e i 



1), 
(9) 

where on , Pi are unobservable phases at low energy. 
The neutrino mass matrix given in Eq. ^ is diago- 
nalized by rh v = W v M v Wl ', and the PMNS matrix 
is determined by W^W^ — £>iVpmnS-D2, where D\ 
and £>2 are diagonal phase matrices similar to those 
in Eq. ©. 



m d [MeV] 
m a [MeV] 
m b [GeV] 

m u [MeV] 
m c [MeV] 
rm [GeV] 

sin^ 3 KM 
sin0f 2 KM 

<5CKM 



1.24 ±0.41 
21.7 ±5.2 
1 nfi+ 014 

l.UD_ 09 

0.55 ±0.25 
210 ±0.21 
82.4 



+30.3 
14.8 



sin 2 0f 3 MNS 
sin 2 0f 2 MNS 

Amli [eV 2 ] 
lAmlxl [eV 2 ] 



0.0351 ±0.0013 
0.0032 ± 0.0005 
0.2243 ± 0.0016 

60° ± 14° 

0.50 ±0.065 
< 0.0155 
0.31 ±0.025 

(7.9 ± 0.3) 10~ 5 



(2-2±°;^) 



10" 



TABLE I: Sample of GUT scale input data used in this work 
(central values and la errors) for Mstjsy = 1 TeV, tan (3 = 
10, and Mqut = 2 x 10 16 GeV (see text for details and 
references). Charged lepton masses are listed in Eq. 11101 . 



B. Input data and \ analysis 

As input data we use quark and lepton masses and 
mixing angles evaluated at the GUT scale, based on 
the RGE analysis of Ref . [35| . As a typical example 
we consider a SUSY scale Msusy = 1 TeV, tan/3 
111 and a GUT scale M GUT = 2x 10 16 GeV. Since 
charged lepton masses are known with an accuracy 
of better than 10 -3 we do not consider them as free 
parameters but fix them to the (GUT scale) central 
values from Ref. 13511 : 



= 0.3585 MeV, 
= 75.67 MeV, 
= 1292.2 MeV. 



(10) 



The remaining data are listed in Tab. D For the 
heavy quark masses m c , mb, nit we adopt the values 
and uncertainties used in Ref. [35[, whereas for the 
light quarks we update the data to the ranges given 
in Ref. m u = 1.5 - 4 MeV, m d = 4 - 8 MeV, 
m s = 80 - 130 MeV (MS masses at 2 GeV). The 
corresponding values at the GUT scale are given in 
Tab.Q] together with the CKM angles and phase [IJ]. 

Concerning the neutrino parameters, we note that 
in the setup under consideration the neutrino mass 
spectrum is generally normal and hierarchical with 
mi < ni2 < ms. In this case the RGE running of the 



PMNS angles [33, [38( is of order 10~ 5 (1 + tan 2 /3) 
10~ 3 , and therefore negligible to a good approxima- 
tion. The running of the neutrino masses is small as 
well (a running effect common to all three neutrino 
masses can be absorbed in the free parameter f v , 



compare Eq. ([6])). In our analysis we use the low 
energy neutrino parameters, as obtained from recent 
global fits to neutrino oscillation data [H, H(| • We 
do not include any constraint on the lightest neu- 
trino mass mi, since the values that are obtained 
are much below the sensitivity of the present data. 

Let us denote the central values and errors of the 
observables by Oi and Ui, where i — 1, . . . , 15 runs 
over all the quantities listed in Tab. HI As described 
above, the predictions of these observables, Pi, de- 
pend on the parameters x a , where a = 1, . . . , 18 runs 
over 3 down-quark masses, 9 real parameters in We, 
fu, fu, \r\, arg(r), arg(^). 1 Then a x 2 -function 
is constructed as 

x 2 (*«)=f:( Pi(Xa l~ 0i Y . (ii) 

i— 1 ^ ' ^ 

The data are fitted by minimizing this function with 
respect to the parameters x a . The minimization 
of the scale factors f u and /„ can be done analyti- 
cally. The remaining 16 dimensional minimization is 
performed with an algorithm based on the so-called 
down-hill simplex method (4lj . 

Let us note that the minimization is technically 
rather challenging. The problem involves para- 
meters which differ by many orders of magnitude, 
and some of the solutions are extremely fine-tuned, 
which leads to very steep valleys in the x 2 landscape. 
Moreover, due to the high dimensionality and the 
non-linearity of the problem there is a large number 
of local minima. In the numerical analysis much ef- 
fort has been devoted to find the absolute minimum, 
involving random methods and dedicated investiga- 
tions to each particular problem under considera- 
tion. However, as it is well known, by numerical 
methods it is difficult to assess with absolute con- 
fidence that an absolute minimum has been found. 
Although each minimum has been carefully tested 
for possible improvements one should always keep 
in mind the possibility that a better solution might 
exist somewhere in the parameter space. 



III. GENERAL FIT OF FERMION MASSES 
AND MIXINGS 

For the discussion of the obtained fits we classify 
the solutions by the values of the parameter |£| which 
controls the relative weight of type-I and type-II see- 
saw mechanism. In Fig. [T] the x 2 minimum is shown 
as a function of this parameter. This means that for 



Since we fix the charged lepton masses to the values given in 
Eq. d 1 1 ) they are neither included in the set of observables 
Oi, nor among the parameters x a . 
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relative contribution of type-I and type-II seesaw 



FIG. 1: x 2 as a function of the parameter |£| which controls 
the relative relevance of type-I and type-II seesaw terms. The 
dashed curve corresponds to the singular (fine tuned) solution 
denoted by mixed' (see text). 



fixed |£|, X 2 ( x a) is minimized with respect to all the 
other parameters x a with a / |£|. In Tab. lU some 
parameter values and the predictions for the data 
are given for four sample points. 

One of the main results of this work is clearly vis- 
ible in Fig. [TJ We find a pronounced minimum for 
|£| ~ 0.36 which corresponds to a mixture of type-I 
and type-II seesaw with a comparable size of both 
terms. Such a "mixed" solution provides an excel- 
lent fit to the data with \ 2 ~ 0.35. From the results 
given in Tab. HTlin the column labeled "mixed" it is 
seen that all observables are fitted within < 0.4er. In 
particular, all the lepton mixing angles and the neu- 
trino mass-squared differences are very close to their 
experimental values. Also CKM CP-violation is de- 
scribed correctly by the value <5ckm = 61° obtained 
in this solution. We conclude that a scenario with 
both seesaw terms of comparable size allows for an 
excellent description of fermion masses and mixings, 
confirming the results of Ref . [25| . Whether this so- 
lution is viable still depends on the detailed study 
of the global vacuum of the given S'O(IO) model. 

The x 2 increases by increasing |£|, and it ap- 
proaches a value of x 2 — 4.3 for |£| > 10, when 
the neutrino mass matrix becomes dominated by the 
type-I seesaw term. Also in the case of complete 
type-I dominance the fit is very good, with most 
observables within < 0.3c, with the sole exception 
of the down-quark mass ma which shows a —1.87a 
deviation from its prediction (see Tab. [Til column 
"type-I"). As we will discuss in more detail later, a 
low value of md is required for a valuable type-I see- 
saw fit. Let us note that all neutrino parameters are 
in excellent agreement with the observations, and in 
particular, the correlation between 0™ NS , 0™ NS , 
and the ratio Ara^/Am^ found in Ref. [25| seems 
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arg(0 

|r| 

arg(r) 
h 

A 



observable 



m d [MeV] 
m s [MeV] 
m 6 [MeV] 

m u [MeV] 
m c [MeV] 
m t [MeV] 

sin^ 3 KM 
sin0f 3 KM 
sin<^ 2 KM 

^CKM [°] 



sm 2 er 3 MNS 

sin 2 0™ NS 
Aml^lO^eV 2 ] 



#PMNS [°] 

« 2 n 



eV* 



type-II 



0.3278 
0.408tt 
16.62 
1.671 x 10~ 10 



pred. pull 



0.7662 
31.33 
1147 

0.5543 
213.1 
78030 

0.0345 
0.00331 
0.2245 
79.35 

0.3586 
0.0145 
0.2829 

7.863 
2.385 
0.279 

-0.70 
1.1 
91 



1.16 
1.85 

0.61 

0.02 
0.17 

-0.29 

-0.43 
0.23 
0.11 
1.38 

-2.17 
0.93 
1.08 

-0.12 
0.50 



14.5 



mixed' 

10" 4 
0.866tt 
1.9977 

I. 849tt 

II. 51 
4.519 x 10" 



pred. pull 



0.4956 -1.82 
22.46 0.15 
1096 0.25 



0.5576 
213.5 
77411 

0.0352 
0.00319 
0.2243 
59.47 

0.5126 
0.0106 
0.3078 

7.894 
2.232 
0.478 

-59 
30 
126 



0.03 
0.18 
-0.34 

0.08 
-0.02 

0.02 
-0.04 

0.19 
0.68 
-0.09 

-0.02 
0.09 



4.1 



mixed 

0.3587 

1.0187T 

0.47896 
0.0013tt 
18.77 
8.732 x 10~ 10 



pred. pull 



1.122 
22.85 
1078 

0.5512 
210.6 
81659 

0.0351 
0.00319 
0.2243 
61.41 

0.5027 
0.0066 
0.3094 

7.898 
2.210 
0.382 

-0.70 
1.8 

-84 



-0.29 
0.22 
0.13 

0.00 
0.03 
-0.05 

0.03 
-0.01 
0.01 
0.10 

0.04 
0.43 
-0.02 

-0.01 
0.03 



0.35 



type-I 

3.59 x 10 6 

1.3187T 

0.3551 
0.0057tt 
19.23 
3.613 x 10" 17 



pred. pull 



0.4719 
19.99 
1029 

0.5538 
213.1 
78117 

0.0349 
0.00323 
0.2243 
61.11 

0.4944 
0.0095 
0.3078 

7.896 
2.223 
0.361 

4.9 
-2.1 
90 



-1.87 

-0.33 
-0.35 

0.02 
0.16 
-0.29 

-0.13 
0.06 
0.01 
0.08 

-0.09 
0.61 
-0.09 

-0.01 
0.06 
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TABLE II: Parameter values and predictions in four example solutions corresponding to different terms dominating the 
neutrino mass matrix: type-I, type-II, or both contributions of comparable size (mixed and mixed'). In the column "pred." 
the predicted values Pi for the observables are given, the column "pull" shows the number of standard deviations from the 
observations, {Pi — Oi)/<Ji, using the data and errors from Tab. [I] Deviations of more than la are highlighted in boldface. The 
final x 2 is the sum of the squares of the numbers in the "pull" column. See the text for comments on the values of the leptonic 
CP phases. 



not to apply here. This solution is not plagued by 
the problem of accomodating the correct <5ckm phase 
found in Ref . [TtJ , and discussed on general grounds 
in Ref. 0. 

For a pure type-II solution with |£| = there 
is more tension in the fit. Although a \ 2 — 14.5 
might be acceptable for 15 data points from a sta- 
tistical point of view, several observables are 1 to 2a 
away from the central value. Among all, one needs 
a large value of the strange-quark mass, namely 
m s ~ 31 MeV that is 1.85cr too large, while the 
PMNS angle 6*23 is too small by 2.2c Furthermore, 
ma, the CKM phase, and the PMNS angle #12 show 
a pull greater than la. These results agree with pre- 
vious analyses of pure type-II solutions [3, HH, EH . 

The dashed line in Fig. [1] corresponds to an inter- 
esting variant of a solution with comparable type- 
I and type-II contributions. Formally this solution 
has a rather small value of |£|, which would signal 
type-II dominance. However, in this case one eigen- 



value of (Md — Mi) is very small, i.e., this matrix 
is close to singular, which implies that its inverse 
has large entries. As a consequence the second term 
in Eq. © turns out to be comparable to the first 
term, in spite of the small value of |£|. Let us de- 
note the type-II term in Eq. ([6]) by M u and the 
type-I term by M l . Then for the solution denoted 
mixed' 'in Fig. Q] and Tab. HTlwe find that the matrix 
entries are of a similar size: 0.2 < \Mjj/MjJ\ < 1.5 
for all i,j — 1, 2, 3. Obviously this solution involves 
a very precise tuning between the Md and Me. matri- 
ces such that the difference becomes close to singu- 
lar. Changing the input values for the down-quark 
masses and the charged lepton mixing matrix We by 
a factor of (1 + 10~ 4 ) destroys in general the fit and 
leads to x 2 values of order 100. 2 Let us add that we 



2 The type-I, type-II, and mixed solutions require tuning of 
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^ 1 c 

^ 15 




m ;/ (M CUT ) [MeV] 



and mixed' solutions given in Tab. [TT] The unshaded region 
corresponds to the lc interval for from Tab. [I] 



obtain mixed' type solutions by extrapolating from 
\(\ = 1(T 4 down to very small values of |£| < 10 9 . 
However, at some point these results are likely to 
become unreliable, since the numerical inversion of 
a nearly singular matrix is limited by the accuracy 
of the algorithm. 



The role of 



charged sector fit is gauged by the need to reproduce 
the tiny electron mass that obeys an approximate 
formula of the form [HI : 

\k'm e \e^ = -|f|e % F d A 4 +e ia2 F c A 6 (12) 
-A 2 A 2 e ia3 X 6 . ^ , , + 0(X 7 ) 



e ia 3 



where m e = m e /m T , A = 1 — p — it] (p, rj being the 
Wolfenstein CKM parameters and A the Cabibbo an- 
gle), while F d = ^/A 4 , F c = ^/A 4 are 0(1) fac- 
tors. The CKM phases cti, /3j are defined in Ref. [HI • 
The parameters f and k! are given by 



-Mr-1\ 



mt 



mt 



fu\r + 3\. (13) 



It is clear that in order to get near the physical value 
m e ~ 2.5 x 10 -4 for \k'\ ~ 0.25 (as suggested by the 
relevant trace identities), the dominant first term on 
the RHS of Eq. (fl~2|) must be either strongly sup- 
pressed (leading to the observed effect of low md 
preference) or cancelled to a large extent by the sub- 
leading ones (often at odds with the CKM phase in 
the first quadrant, c.f. [0],[24|]). Therefore, for low 
values of the available portion of the parameter 
space is larger and allows for a better global fit of the 
remaining physical parameters. It is perhaps worth 
mentioning that there is another pattern in our data 
(though by far much weaker than the low mj pref- 
erence) that can be justified on the same grounds, 
namely a generic drift towards lower m* and higher 
m c regions. One verifies that in such a case the sub- 
dominant term proportional to F c gets larger and 
allows for a better "screening" of the dominant first 
term on the RHS of formula ([121). 



The down-quark mass plays a relevant role in ob- 
taining a good fit to the data. As it is visible from 
Tab. ITTlthe type-II, mixed', and type-I solutions re- 
quire values of which are about 1.2c, 1.8c, and 
1.9a below the preferred value given in Tab. [T] In 
Fig. [2] we illustrate how the fits become worse when 
rrid is increased. Especially the very good fits of 
the type-I and mixed' solutions with \ 2 ~ 4 are 
strongly affected by increasing ma, and values of 
m d ~ 1.5 MeV at M GUT lead to \ 2 - 15 and 20, re- 
spectively. Also for the pure type-II case the fit soon 
becomes unacceptable if md is increased. Only the 
mixed solution has enough freedom to accomodate 
larger values of the down-quark mass, and in this 
case also a very good fit is possible for ~ 2 MeV. 

As a matter of fact, it is quite easy to under- 
stand these results on an analytical basis. As it was 
pointed out in Rcfs. [lj] and [2J] the quality of the 



the parameters with a typical accuracy of better than 10 3 
(the tuning is slightly less severe for the type-II case). 



B. Predictions for the neutrino parameters 

In this section we discuss in some detail the pre- 
dictions of the studied setup for the neutrino sector. 
Our main results are summarized in Fig.[3l where we 
show how the x 2 changes if 6»p mns , 6>f 3 MNS , and the 
mass of the lightest neutrino mi are varied. Techni- 
cally this analysis is performed in the following way: 
to test the variation of an observable Ok the term 
with i = fc is removed from the % 2 given in Eq. (|lip 
(this ensures that the only constraint on the observ- 
able comes from the mass sum rules and not from 
the input data). Then, in order to test a certain 
value O* for the observable Ok a term is added to 
the x 2 with a very small error of 1%, to confine the 
fit: (P k (x a )~O*) 2 /(0.01O*) 2 . After the minimiza- 
tion this term is removed, and the x 2 is evaluated at 
the point obtained in the minimization. 

Fig. [3] (left) shows the constraint on the PMNS 
mixing angle #2 3 MNS - One can see that for the type- 
I, mixed, and mixed' solutions there is no definite 
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FIG. 3: x 2 as a function of sin 2 0f 3 MNS (left), sin 2 6»f 3 MNS (middle) and of R = m^Am^ (right) for the type-I, type-II, 
mixed, and mixed' solutions given in Tab. |n] The shaded regions are excluded at 2a according to the data given in Tab. U 
The dotted vertical line shows roughly the sensitivity to $P g MNS of neutrino oscillation experiments within a timescale of 10 
years [43l |. 



prediction for this angle and very good fits are pos- 
sible with values of #2 3 MNS in the whole range al- 
lowed by the data. In contrast a scenario with pure 
type-II seesaw shows a clear preference for small val- 
ues of sin 2 #23 MNS - I n particular, maximal mixing 
= 0.5 is disfavoured with respect to the 



sin 2 0P 3 MNS 



best fit value sin" 2 02g MNS 



0.36 with A X A 



11. 



Hence, the pure type-II seesaw model predicts siz- 
able deviations from maximal mixing within the 
reach of upcoming neutrino oscillation experiments, 
see, e.g., Refs. lio. |42|. This result is in agreement 
with Refs. 0, |25| . Let us stress, however, that de- 
viations from maximal mixing are not a general pre- 
diction of the S'O(IO) model under consideration; it 
holds only for the pure type-II case. 

Concerning the mixing angle fl™ 1 ^, for all solu- 
tions the best fit point predicts values close to the 
present upper bound and clearly within the reach of 
upcoming neutrino oscillation experiments H3, 5H . 
Also in this case the pure type-II solution gives the 
most stringent prediction. However, if the fit is 
stretched to some degree, the type-I and mixed so- 
lutions allow also for smaller values of 0™ that 
might be difficult to detect in the next round of ex- 
periments. For instance, in the mixed solution case a 
fit with x 2 ~ 7 is possible for sin 2 6>f 3 MNS = 2 x 10~ 3 . 
An interesting feature of the mixed' solution is that 
even for very tiny values of sin 2 #f 3 MNS < 10~ 4 a 
reasonable fit can be obtained with x 2 ~ 12.6. The 
main contributions to the x 2 i n this case are devia- 
tions of —1.7cr for md, —1.2a for m s , —2. Oct for m t , 
and +1.4ct for sin 2 #23 MNS . The existence of this so- 
lution shows that if no signal of 6*f 3 MNS is detected 
by the upcoming experiments it might be still pos- 
sible to construct models with viable predictions for 



fermion masses and mixings, although the amount 
of fine tuning will increase. 

Another interesting information the setup under 
consideration provides, concerns the shape of the 
neutrino mass spectrum. The most solid predic- 
tion is the normal mass ordering, which means that 
mi < nri2 < m 3 , i.e., Am§ x > 0; no viable solution 
has been found for inverted ordering (Am 31 < 0) 
which is equally allowed by present oscillation data. 
The predictions for the absolute neutrino mass scale 
are given in terms of the ratio of the lightest neu- 
trino mass m% to the square root of the "solar" mass- 
squared difference: 



R 



mi 



(14) 



The best fit values for this ratio given in Tab. [Hi for 
the four example solutions are in the range 0.28 < 
R < 0.48. These values show that there is only a 
modest hierarchy in the neutrino masses. For exam- 
ple, a value of R — 0.3 implies that m,2 — 3.5 mi, 
i.e., mi and mi are of the same order of magnitude. 
From the plot in the right panel of Fig. [3] one can 
infer the allowed ranges for the ratio R. We find 
reasonable fits for values in the range 



0.2 <R<2. 



(15) 



where the upper bound implies mi ~ 1.1 mi. For 
the pure type-II solution the ratio R is stronger con- 
strained to values around the best fit point of 0.28, 
whereas the constraint is weakest for the mixed so- 
lution. Note that a quasi-degenerate neutrino spec- 
trum with mi ~ m .2 ~ m.3 would correspond to 
R > yj AmyAm^ ~ 5.3 which is clearly excluded 
within the 50(10) framework under consideration. 
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The range for R given in Eq. (fT"5|) implies the fol- 
lowing intervals for mj and the sum of the neutrino 

masses E = mi + m-2 + m^: 



1.8 x 1(T 3 eV < mi < 1.8 x 10~ 2 eV . 
0.058 eV< E< 0.088 eV. 



(16) 



For comparison we note that for mi ~ one has 
S = 0.056 eV, which shows that the smallest values 
of mi, for which reasonable fits are found, are close 
to the case where the contribution of mi to E can 
be neglected. 3 

For any given set of input parameters the val- 
ues of the phases in the PMNS matrix are deter- 
mined as well. We have investigated in the four cases 
of Tab. Hi] the predictions for the Dirac CP phase 
<5pmnSj as well as the two Majorana phases a,\ and 
«2j defined in analogy with Eq. ([9]) in the leptonic 
sector. In all cases we find a correlation among the 
phases, such that 



Jpmns 



-1(X\ ~ — 2(I 2 + 7T. 



(17) 



For the mixed, mixed', and type-I solutions there is 
no definite prediction for the values of the phases, 
and, given the correlations above, fits of comparable 
quality are found for all values of the phases in the 
[0, 2tt] range. Only in the case of pure type-II see- 
saw the fit prefers values of Jpmns ~ 0, ot\ ~ (ir), 
ct2 ~ ±tt/2, in agreement with the results of Ref. [2J] 
(the Majorana phases 0i,2 there reported are defined 
according to Ref. [37j as —2a\^)- Moving (Jpmns 
from the preferred value to ±60° increases the \ 2 by 
10 units. 



126// responsible for the type-II seesaw is induced 
via 210// couplings as well. The model is minimal in 
the number of parameters, 26 altogether (soft SUSY 
breaking aside). The same number is found in the 
minimal supersymmctric standard model (MSSM) 
with right-handed neutrinos and it is much less 
than the number of parameters in the correspond- 
ing SUSY 5t/(5) GUT. 

Since none of the 126// nor 210// Higgs multiplets 
can couple to the 5O(10) matter bilinear 16 m (g) 16 m 
at the renormalizable level, the Yukawa superpoten- 
tial reads 

W Y = 16 m (Yi 10// + Yi 26 T26//)16 m , (18) 
while the 5O(10) Higgs sector is described by 



W H = 210// 2 



4! 



210 



Mi 



26 



H 



5! 



126// 126// 



~126 H 210// 126// + MiqIO// 2 
5! 



210// 10// {a 126// + a 126//) . (19) 



The minimization of the scalar potential has been 
analyzed in great detail in Refs. [HI, H^, [3l|. The 
need of a careful study of the GUT scale poten- 
tial in order to perform consistent predictions of the 
fermion mass textures is emphasized in these papers. 
The light MSSM- like Higgs doublets entering the low 
energy Yukawa potential are in general a superposi- 
tion of corresponding doublet components of all the 
Higgs multiplets. The set of VEVs which describes 
the vacuum of the model is governed by one complex 
parameter x fl5j ] : 



IV. THE MINIMAL SUSY SO(10) GUT 



Having analyzed the implications on the fermion 
mass fit of the 10// plus 126// 5*0(10) Yukawa sector 
we focus now on the study of the so called minimal 
renormalizable SUSY 5O(10) GUT [l|. The model 
is characterized (in addition to 10// and 126//) by 
the presence of the 126// and 210// representations 
in the Higgs sector. The 126// is needed to pre- 
serve the GUT-scale D-flatness while 210// plays the 
dual role of triggering the spontaneous 5*0(10) gauge 
symmetry breaking and provides the necessary mix- 
ing among the 10// and 126// weak doublet com- 
ponents needed to achieve (after electroweak sym- 
metry breaking) a realistic fermion spectrum. The 
tiny VEV of the left-handed triplet component of 



3 The fact that the fit always gives a hierarchical neutrino 
mass spectrum justifies a posteriori the neglect of the run- 
ning of the neutrino parameters [37l |38| . 



\1> 1: 1)210 — 



(15,1,1) 



M210 x(l - 5x 2 ) 
A (l-x) 2 ' 
M210 \l-2x-x 2 ) 



A 



(1-x) 



(20) 



(15, 1, 3) 210 — - x , 

/— \ x 2M 9 2 ln x(l - 3x)(l + x 2 ) 

<io, 1, 3> 126 (10, 1, 3) w = —m [ n _^ a > 



{l-xf 



Avoiding D-tcrm SUSY breaking requires 
<T0,1,3) 126 = (lO,!,^ = ««■ In Fig. 13 we 
depict the relevant 5O(10) breaking patterns 
together with the related VEVs. 

The various mass parameters in the GUT super- 
potential can be written as functions of x as well. 
For instance 

f 3 " Ux + 15x2 - 8x3 
M 126 = A/210 I • (21) 

The mass parameter M 10 is determined by the min- 
imal fine-tuning condition (l3l . [l5[ which drives the 
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FIG. 4: Symmetry breaking patterns and related VEVs in the minimal SUSY 50(10). We do not display breaking chains 
involving intermediate SU(5) symmetries, at odds with the proton decay constraints. Pati-Salam notation is used. 



mass of two Higgs doublets down to the weak scale. 
The relevant formula reads 



Mio = M 2 



Pw 



2-q\ (x - l)p 3 P5 ' 



(22) 



where p n are polynomials of x, see Refs. [la, |29|. 
The SU(2)l triplet mass, relevant for the type-II 
seesaw, can be written as 

V x(Ax 2 - 3x + 1) .„„. 
M T = M 210 1 — 23 

a (x — \y 

while the tiny induced SU(2)l triplet VEV is given 
by [H,[3l| 



v L = (10,3,1)™ = 



(ay™ + V6f]vj 26 )v 



210 



(24) 



The electroweak VEV components v u 's are deter- 
mined in terms of the electroweak scale and the x 
parameter, c.f. Refs. [T5L [29t| - 

In the first part of the paper we learned that the 
best fit of the present data on fermion masses and 
mixing is obtained by a mixed contribution of type-I 
and type-II seesaws. The model independent analy- 
sis based only on the form of the Yukawa potential 
in Eq. (| 18(1 assumed the vacuum state to provide 
the required set of parameters, in particular the cor- 
rect neutrino mass scale. In the following we exam- 
ine thoroughly the question whether such a vacuum 
exists in the minimal renormalizable SUSY SO(10) 
framework. The crucial issues are gauge coupling 
unification and the present limits on proton decay. 
Both instances in a low energy SUSY setting priv- 
ilege a great desert scenario up to 10 16 GeV. This 
requirement in turn tends to largely suppress the 
neutrino mass scale as obtained from type-I and/or 
type-II seesaw. 



A. Matter fermion fit in the minimal 
renormalizable SUSY 5O(10) model 

For our numerical analysis we adopt the parame- 
terization and the phase convention of Ref. [29j . Af- 
ter minimal fine-tuning the Higgs potential can be 
described in terms of 8 real parameters: 

M 2W , Re(x), lm(x), a, a, |A|, \r)\, 

= arg(A) = -arg(?7) . (25) 

For our purpose (fermion masses and symmetry 
breaking) without loss of generality we can set <fi = 0. 
As emphasized in Refs. [HlHIj], the parameter x pro- 
vides a systematic and effective way of describing the 
5O(10) symmetry breaking patterns and the related 
fermion mass scales. 

In the minimal SUSY 5O(10) framework under 
consideration the VEVs appearing in the fermion 
sum rules Eq. (fTJ) are no longer free parameters 
(as assumed in the numerical analysis of Sec. IIII[) . 
but are functions of the Higgs potential parameters, 
c.f. Eq. (|23|) . Using the results and notation of 
Ref. [29j one obtains the following expressions for 
the VEV combinations appearing in the fermion fit 
as defined in Eqs. and ([7]): 



f u = ^tan/3- 



fu 



J\I 2 



■ tan {3 sin j3 a 




16 fn(x) " 



(26) 



Here N u , are functions of the parameters in Eq . 
([25]) given explicitly in the appendix of Ref. [29l |. 
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FIG. 5: x 2 contours of the fermion mass fit in the x plane for the sample data in Table I. All scalar quartic couplings are set 
to one. The neutrino mass scale factor /„ is treated as a free parameter. The black contour curves enclose the regions where 
f v (x)/f v {?it) > 0.005, where f v (x) is calculated from Eq. II26I I and /„(fit) is the value obtained from the fermion fit. The cross 
symbols corresponds to points where either fj or fjj has a singular behaviour. 



v = 174 GeV is the electroweak scale, while 



2 P2P5 , (x- l)(4x- l>3g|g 
III = 7, 



P3CT 



where 



'2x(l - 3x)(l + x 2 ) 



(28) 



and p n , Qn are polynomials of x of order n (see Tab. I 
in Ref. |5|). 4 

We now address the question of whether a real- 
istic fermion fit is possible, given the constraints 
from the Higgs sector encoded in Eqs. (f2"o| which 
together with Eqs. ©, Q, and © lead to a highly- 
constrained system of relations. First we note that 
not much freedom is left in adjusting the parameters 
a, a, r], A. This follows from the functional form of 
the VEVs in Eqs. (|2"0"|) . of the potential mass pa- 
rameters, and of f v ,N u ,Nd, as well as from the 
requirement of perturbativity of the Higgs poten- 
tial. Hence, in what follows we set a ~ a ~ r\ ~ 
A ~ 0(1) 29\. Similar considerations suggest also 
|x| ~ 0(1), and therefore we first restrict our search 
to the range 0.1 < |x| < 10 (deferring a discussion 
on the small and large x regimes to the end of the 
section). 



In order to identify candidate solutions we per- 
form a scan in the complex x plane by first fixing 
a = ot = i] = \ = 1. The x dependence of r and 
£ is taken into account according to Eq. (|26p. leav- 
ing however the neutrino mass normalization /„ as 
a free parameter to be determined from the fit. Also 
f u becomes a free parameter, if tan/3 is allowed to 
vary. 5 

In Fig. [5] we show the % 2 contours of the fermion 
mass fit in the complex x plane, where darker ar- 
eas represent better fits of the fermion masses and 
mixings according to the % 2 values given. As it is 
visible from the figure we find large regions in the x 
plane where a very good fit with x 2 < 10 can be ob- 
tained. The crucial question is whether within these 
regions it is possible to obtain also a consistent value 
of f v . From Eq. ([6]) one finds that /„ or f v \£\ has 
to be of order \J Am^/mj ~ 5x 10~ n to provide a 
neutrino mass scale required by oscillation data, in 
agreement with the values reported in Tab. [ill In 
contrast, from Eq. (|26p one finds a natural size of 



fu 



utan/3 
M 2W 



5 x 10 



-13 ( tan/3 
V 55 



(29) 



for M 2 io = 2 x 10 16 GeV. This estimate illustrates 



4 Notice that our parameter £ controls the relative weight of 
type-I and II seesaw, and is different from the £(x) used in 
Ref. [H. 



5 In the following we will consider tan f) as a free parameter 
within the range 10 < tan f) < 55. We may still keep the 
data in Table [I] as reference GUT values for the fermion 
parameters, since changing tan fi in the indicated range has 
a small impact on the running [35| , 
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observable 



m d [MeV] 
m s [MeV] 
m b [MeV] 

m u [MeV] 
m c [MeV] 
m t [MeV] 

sin^ 3 KM 
sin0? 3 KM 
sin0? 2 KM 

<5cKM [°] 

sin 2 ^ 3 MNS 

sinX 3 MNS 
sin 2 0f 2 MNS 

Aml^lO-^V 2 ] 
AmlillO-^V 2 ] 



mil 



parameter 


best fit 


A, a, a fixed 


\x\ 


1.0000002207 


1.00000009837 


arg(x)/7T 


0.4999999781 


0.4999999815 


?7 


0.3191 


0.3191 


A 


1.5642 


1 


a 


2.9455 


1 


a 


4.0332 


1 


M 2 io [GeV] 


2.000 x 10 16 


2.000 x 10 16 


tan/3 


41.15 


45.88 


Kl 


7.476 x 10 4 


1.513 x 10 5 


arg(C)/ 71 " 


-0.60492 


-0.53116 


|r| 


1.955 


1.955 


arg(r) /it 


-0.2568 


-0.2568 


fu 


12.28 


11.47 


u 


7.077 x 10~ 17 


3.745 x 10" 17 


(1, 1, 1) [GeV] 


3.836 x 10 16 


6.000 x 10 1B 


(15, 1, 1) [GeV] 


2.557 x 10 16 


4.000 x 10 1B 


(15, 1, 3} [GeV] 


1.279 x 10 16 


2.000 x 10 1B 


(10, 1, 3) [GeV] 


3.423 x 10 13 


3.009 x 10 13 


M PG [GeV] 


4.135 x 10 10 


2.043 x 10 10 


M T [GeV] 


8.655 x 10 15 


1.354 x 10 16 



pred. pull 



0.3435 
27.91 
1065 

0.5596 
212.8 
82090 

0.0351 
0.00325 
0.2244 
64.32 

0.4893 
0.0133 
0.3007 

7.875 
2.199 
0.3351 



2.2 

1.2 

0.038 

0.038 
0.15 
-0.21 

-0.039 
0.10 
0.033 
0.31 

-0.16 
0.86 
-0.37 

-0.082 
-0.002 

7.25 



pred. pull 



0.3373 
26.36 
1112 

0.5639 
213.1 
79550 



2.2 

0.90 
0.37 

0.056 
0.16 
-0.19 



0.0350 -0.065 

0.00329 0.18 

0.2244 0.057 

70.10 0.72 

0.5002 0.0035 
0.01085 0.70 
0.2953 -0.59 

7.860 -0.13 
2.274 0.20 

0.3289 

7.30 



The only way left to obtain the needed neutrino 
mass scale is to consider special points in the x plane, 
where either or |//(a:)| becomes large be- 

cause of roots in the denominators of Eq. |27p . These 
singularities are marked in Fig. [5] by the crosses. 
The closed black contours correspond to the regions 
where f v (x)/f u (&t) > 0.005. We find that only close 
to the singularities /„ can be large enough to pro- 
vide the required neutrino mass scale. On the other 
hand, it appears that most of the singularities drop 
into regions where the fermion mass fit is poor. We 
carefully checked numerically the regions close to 
the relevant singularities and came to the conclu- 
sion that only in the near neighborhood of x = ±i a 
realistic fit of the matter sector is possible. Helpful 
discussions of specific points and limiting cases in the 
x complex plane based on analytical considerations 
can be found in Refs. [l5|, H^] . 

In Tab. lIIII we report the best fit solution we found. 
The absolute value and phase of x are tuned at the 
level of 10 -7 in the close neighborhood of i. The 
"symmetric" solution x = — i exhibits to quite simi- 
lar features. In the "best fit" column we have opti- 
mized the parameters a, a, rj, A, whereas in the third 
column we have taken a = a = X = 1. As expected, 
changing these parameters has very little impact on 
the fit, as it is evident from the comparison of the 
second and third column. Requiring in addition 
rj = 1 gives a very similar result. The neutrino mass 
matrix is dominated by type-I seesaw since a = 
for x — ±i, which implies // — » 00 and fu = 0, and 
hence £ — > 00, see Eqs. (127)) and (|26p. From the data 
in the table one finds /„|£| ~ 10 -11 , which is the 
correct value to provide the required neutrino mass 
scale. The table shows that a very good fit of all 
quark and lepton mass and mixing parameters (in- 
cluding the neutrino sector) is obtained, with just a 
2.2(7 pull for TTid that is acceptable in view of the sys- 
tematic uncertainities entering the light quark mass 
determinations . 



TABLE III: Parameter values and outcomes for the best fit 
solution, x ~ i. The reference mass parameter M210 is set at 
a typical MSSM GUT scale. The charged lepton masses are 
given in Eq. lit Oi l. Deviations above 2<r in the fermion mass 
data are highlighted in boldface. The final x 2 is the sum of 
the squares of the numbers in the "pull" column. 



that in the framework under consideration a gap 
of about 2 orders of magnitude exists between the 
generic prediction of the neutrino masses and the 
scale required by the data. We have verified this 
expectation numerically and for most part of the 
complex x plane shown in Fig. [5] the value of /„ (x) 
according to Eq. (|26p is 2 to 4 orders of magnitude 
smaller than the value of f u required by the fit (de- 
noted by /„(fit)). 



B. Gauge couplings unification and proton 
decay for x ~ ±i 

What remains to be checked is whether the best fit 
solution satisfies unification and proton decay con- 
straints. For x — ±i the 5*0(10) gauge symmetry 
is broken 6 to G3211 = SU{3) C X SU(3) L xU(l) R x 
U{1) B -l bythe VEV of (15, 1, 3) 2 io (compare Fig.g) 
and Refs. [H, H§|). The calculation of the particle 



6 We neglect the tiny hierarchy among the VEVs of 
(1, f , 1)2101 (15, 1, 1)2101 and (15, 1, 3)210 a $ all of them are 
confined for x = ±i within an interval of less than one order 
of magnitude, c.f. Eq. (20) and Tab.lTTll 
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M PS 




FIG. 6: Running gauge coupling constants in the setup of 
Table ITTll We have M SU SY = 1 TeV, M PG = Af(8,3,o)> 
M3211 = (10,1,3), M 322 i = (15,1,1), M PS = (1,1,1), and 
Mqut = M210 = 2 x 10 16 GeV. 



gauge couplings unification. The minimal renormal- 
izable SUSY 5*0(10) scenario seems therefore fail- 
ing to provide a realistic description of the low en- 
ergy world. In general, if right-handed neutrinos 
trasform non-trivially under the GUT gauge sym- 
metry, the minimum neutrino mass scale required 
by oscillation data implies via type I (or type II) 
seesaw the presence of states at intermediate mass 
scales [29], [3(J, that spoil the successful unification 
within the MSSM. 

In passing let us stress that the example discussed 
above shows quite clearly that the naive argument 
of setting all relevant particle states at the scale of 
the symmetry breaking step may be far from cor- 
rect. One should always be aware that much lighter 
states may appear in the spectrum as a consequence 
of the spontaneous breaking of accidental (would-be) 
symmetries of the scalar potential 



spectrum shows that in this limit a set of unmixed 
states that transform as (8,3,0,0) under G3211 (or 
(8,3,0) under the SM gauge group) remains light. 
In the exact limit x = ±2 these Goldstone states 
are massless because of the spontaneous breaking of 
accidental global symmetries of the Higgs potential. 
Upon breaking the G3211 gauge symmetry to the SM 
via the VEVs 7 (10, 1, 3)^ = <10, 1, 3> 12fl = v R (= 
-M3211) they acquire a mass Mpq of the order of 
M§ 211 /M PS , where M PS = (1,1,1) 210 is the Pati- 
Salam scale. According to Ref. [15| the expression 
for the Pscudo-Goldstone mass is 



M PG = -4M; 



210 ■ 



(2x- l)(x 2 + 1) 
(x- l) 2 



(30) 



which is exactly zero for x = ±4. In our best fit case 
reported in Tab.|ffl]we find M 32 n ~ 10 13 GeV, while 
Mp S ~ 10 16 GeV. One therefore expects M PG ~ 
10 10 GeV, as we consistently find (see table). 

In spite of the fact that these states are much 
lighter than the GUT scale they do not affect proton 
decay because of their zero B — L charge. However, 
they do affect heavily the running of the gauge cou- 
plings (they transform as adjoint representations of 
SU(3) C and SU(2)l). Fig.[H]shows how dramatically 
the running is affected. The unification of the gauge 
couplings is completely spoiled, with the coupling 
constant of SU(2)l diverging below the GUT scale. 
Hence, the only solution found providing a realistic 
fit to all fermion masses and mixing parameters has 
to be discarded because of the dramatic failure of the 



7 More precisely, the quantum numbers of the states 
responsible for this symmetry breaking step are: 
(l,l,±l,q=2) 3211 £ (1,1,3,T2) 3221 £ (10, 1, 3)^ © h.c 



C. Large/small \x\ regime 



Before conclusively rejecting the minimal renor- 
malizable SUSY 5*0(10) scenario we must make sure 
that the domain 0.1 < \x\ < 10 considered in the pre- 
vious analysis covers the physically relevant region. 
Indeed, the shape of the x 2 contours in Fig. [5] sug- 
gests that all the matter fermion data might be well 
reproduced in areas outside the considered region. 
Let us give simple semianalytic arguments against 
this option. 

By observing that in the large x regime all 210# 
VEVs in Eq. (|20|) scale as M 2 io \x\/X and proton 
decay forces us to maintain most of the GUT states 
above Mqut — 10 16 GeV, we must require 



Mgio N 
A 



10 16 GeV 



(31) 



In the large x limit (and for large tan/3) from 
Eqs. PSI) . (|2T|) and the expressions for N u , N d given 
in Ref. [29( one finds for the neutrino mass scales 

MC\ oc - i ; tan/3, f v oc , (32) 

M 2 w\x\ M 2 ioF| \ x \ 

in the type-I and type-II seesaw cases respectively. 
By requiring M 2 io|a;| — const (see Eq. (13TT0 these re- 
lations show that also for large x the neutrino mass 
scale cannot be efficiently enhanced over its natu- 
ral value i> 2 /-Mgut, and accordingly our numerical 
calculation yields values of /„ of about four orders 
of magnitude too small in the region \x\ > 10. In 
passing let us remark that in the |x| S> 1 regime 
there appear a number of states with masses around 
M210 < Mqut that again spoil gauge unification, 
albeit not affecting proton decay (see Ref. [Hj]). 
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In the small \x\ regime one finds for both seesaw 
types 

/„,/»,|£|oc— Vr^^VR- (33) 

M 2W \x\ 

Once again taking into account Eqs. (|2QI3ip and the 
fact that scalar couplings are bound to vary in a 
narrow range by perturbativity on one side and mass 
scale splittings on the other, no substantial enhance- 
ment of the neutrino mass can be expected. Our x 2 
fit does confirm numerically these expectations. 

We therefore conclusively show, independently 
confirming the conclusions of Ref . [3(| , that the min- 
imal renormalizable SUSY 50(10) setup in spite of 
noteworthy and impressive features in reproducing 
the observed flavor textures fails in reproducing the 
neutrino mass scale. Generally, in a renormalizable 
framework, type I and type II seesaws call for inter- 
mediate scales that spoil (SUSY) gauge unification. 



V. CONCLUSIONS 

In the present paper we studied in great detail 
the fit of quark and lepton masses and mixing data 
within a class of supersymmetric 5*0(10) grand uni- 
fied models with a minimal renormalizable Yukawa 
sector based on one 10h and one 126// Higgs rep- 
resentation. A systematic optimization of the rele- 
vant x 2 function has been performed, complemen- 
tary to random parameter searches applied in previ- 
ous studies. We have shown that for a comparable 
size of the type-I and type-II seesaw terms, an excel- 
lent fit to the fermion data can be obtained, where 
all observables (including the CKM CP phase) are 
fitted within less than 0.4 standard deviations. So- 
lutions based on pure type-I and type-II seesaw have 
been discussed as well, and we have identified a new 
class of possible solutions corresponding to a singu- 
lar behaviour of the type-II Majorana mass matrix. 
The corresponding predictions for the neutrino pa- 
rameters have been investigated in detail. 

In the second part of the paper these general re- 
sults were confronted with the additional restric- 
tions emerging in the minimal renormalizable SUSY 
5O(10) scenario from the model vacuum as well as 
from proton decay and unification constraints. We 
identified a very limited (and fine-tuned) area in 
the parameter space where, given the constraints 
from the minimization of the Higgs potential, a fit 
to the quark and lepton masses and mixing para- 
meters is possible. However, all solutions provid- 
ing a realistic fit to the matter fermion data had to 
be discarded because of the failure in reproducing 
the strenght of the low energy gauge couplings. In 



GUTs that embed non trivially right-handed neu- 
trinos the absolute neutrino mass scale required by 
oscillation data generally implies via type I and/or 
type II seesaws the presence of states at intermedi- 
ate mass scales, that are likely to spoil the successful 
unification achieved in the MSSM. Our analysis did 
in fact exclude, in the case of the minimal SUSY 
50(10) model, the existence of (fine tuned) excep- 
tional solutions. 

The minimal renormalizable SUSY 50(10) sce- 
nario based on one 10//, one 126//, one 126//, and 
one 210// Higgs representations is conclusively re- 
jected. Our discussion suggests that either one at- 
tempts to enlarge the scalar (Yukawa) potential in 
order to soften some of the parameter correlations 
on the vacuum, or non-renormalizable terms are in- 
cluded providing effectively the needed depletion of 
the seesaw scale. Alternatively, one is let to con- 
sider the intriguing option of a non-supersymmetric 
GUT scenario, where intermediate scales, preferred 
by neutrino data, are needed anyway by gauge cou- 
pling unification. Beauty and loss do strike our deep- 
est chords. As the poet says: 
Is that time dead? lo! with a little rod 
I did but touch the honey of romance — 
And must I lose a soul's inheritance? 
(Oscar Wilde, Helas, 1881). 
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Note added 

The arXiv version 4 of this paper corrects three 
entries in Table III, namely the value of r\ in the third 
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two values of f v that were off by a factor two. We 
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for carefully checking our results and spotting the 
typos. 
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